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I.  INTRODUCTION 


Since  1972,  Applied  Research  Laboratories,  The  University  of  Texas 
at  Austin  (ARL : UT ) ,  has  been  heavily  involved  in  research  and 
development  in  low  frequency  ocean  acoustics  for  the  U.S.  Navy.  Most 
of  this  work  has  been  exploratory  and  advanced  development.  It  has 
included  development  of  numerical  modeling  for  both  research  and 
applications;  i.e.,  planning,  collection,  and  analysis  of  planned  and 
existing  acoustic  systems.  One  aspect  of  this  ARL:UT  effort  is  that  it 
exposes  questions  which  are  appropriate  to  basic  research,  and  for  which 
answers  will  find  ready  applications. 

Part  of  the  overall  propagation  studies  at  ARL : UT  is  a  basic 
research  program  funded  by  Office  of  Naval  Research  (Code  425)  devoted  to 
measurement  and  modeling  of  acoustical  properties  of  ocean  sediments. 
The  research  program  originated  as  an  effort  to  measure  acoustical 
parameters  of  sediments  in  situ  during  coring  operations.  An  instrument 
called  a  profilometer,  described  in  many  of  the  references  contained  in 
Appendix  B,  was  developed  to  measure  compressional  wave  velocity,  shear 
wave  velocity,  and  impedance  profiles.  Use  of  the  profilometer  requires 
only  minor  changes  to  the  cutter  head  of  a  geophysical  corer;  it  has 
been  used  on  several  cruises. 

During  the  development  of  the  profilometer,  it  was  necessary  to 
improve  the  sensitivity  and  size  of  acoustic  transducers.  In 
particular,  this  work  provided  a  significant  improvement  in  shear  wave 
transducer  sensitivity;  a  patent  was  issued  for  this  design.  These 
transducer  improvements  proved  to  be  just  as  useful  for  basic  laboratory 
measurements  as  for  field  work.  Thus,  part  of  this  program  began  to 
include  laboratory  work  with  carefully  controlled  sediments,  and 
theoretical  analytical  work  as  part  of  the  interpretation  of  those 
measurements. 
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During  1981,  the  program  began  to  discontinue  development  of  the  in 
situ  profilometer  and  to  concentrate  on  theoretical  analysis  and  labora¬ 
tory  measurements.  To  bolster  the  theoretical  effort,  faculty  and 
students  of  the  Department  of  Aerospace  Engineering  and  Engineering 
Mechanics  of  The  University  of  Texas  at  Austin,  who  specialize  in  the 
study  of  mechanical  properties  of  composite  materials,  are  participating 
in  this  program  on  a  regular  basis. 

During  1980,  laboratory  measurements  were  begun  on  artificial 
(glass  bead)  sediments  with  different  mixtures  of  glycerol  and  water  as 
a  pore  fluid.  During  1981,  these  results  were  analyzed  and  prepared  for 
publication  (Appendix  A  of  this  report  contains  a  preprint  of  this  pub¬ 
lication).  Then,  to  clarify  some  points  which  arose  during  that 
analysis,  measurement  techniques  were  modified  and  the  experiment  was 
repeated.  Those  results  are  included  in  Section  III  of  this  report. 

A  bibliography  of  publications  under  the  sediment  acoustics  program 
is  included  as  Appendix  8.  Since  the  program  was  started,  15  technical 
reports  have  been  published,  14  papers  have  been  presented  at  technical 
meetings,  eight  papers  have  been  published  in  or  presented  to  scientific 
journals,  three  papers  have  been  included  in  books,  and  two  invention 
disclosures  have  been  submitted  for  patent.  Of  these,  in  1981  one  tech¬ 
nical  report  was  issued,  one  publication  submitted,  and  one  technical 
presentation  given. 


II.  SEDIMENT  ACOUSTICS  ISSUES  IN  LOW  FREQUENCY 
SOUND  PROPAGATION  IN  THE  OCEAN 


A.  Low  Frequency  Propagation  and  the  Seabottom 

As  mentioned  in  the  introduction,  ARL:UT  is  involved  in  several 
areas  of  low  frequency  ocean  acoustics.  This  exposure  has  pointed  out 
several  issues  which  we  feel  will  be  supported  at  a  basic  level  by  the 
research  under  this  program.  In  this  section,  we  first  briefly  identify 
some  of  these  problems.  Then  we  discuss  some  of  the  issues  related  to 
theoretical  analysis  of  the  acoustical  properties  of  sediments. 

Studies  of  low  frequency  (25-1000  Hz)  ocean  acoustics  are  usually 
divided  between  deep  ocean  regions  (depths  of  2000  m  or  more)  and 
shallow  water  regions  (depths  of  2000  m  and  less).  Usually,  depths  of 
500-2000  m  are  continental  slopes  and  are  thus  localized;  for  that 
reason,  500  m  is  a  practical  maximum  depth  for  shallow  water.  A  general 
study*  of  the  bottom  acoustics  issues  for  exploratory  development  work 
has  recently  been  prepared  by  Naval  Ocean  Research  and  Development 
Activity  (NORDA)  Code  530;  it  goes  into  some  detail  about  acoustical 
mechanisms  and  their  relationship  to  Navy  acoustical  systems. 

For  deep  water  propagation,  the  quantity  of  most  interest  is  the 
loss  of  energy  per  bottom  interaction,  i.e.,  per  refraction  in,  or 
reflection  from,  the  bottom.  That  quantity  has  in  the  past  been 
characterized  as  "bottom  loss."  However,  the  more  sophisticated  current 
applications  model  the  bottom  as  part  of  the  propagation  medium,  and 
thus  require  an  understanding  of  bottom  acoustics.  Deep  ocean 
sediments,  with  small  grain  sizes  and  soft  compositions,  are  readily 
penetrated  by  sound.  Examining  typical  velocity  gradient  and 
attenuation  values,  one  can  deduce  that  the  properties  of  sediments  at 
200-500  m  depth  are  important  at  50  Hz.  At  200  Hz,  depths  greater  than 
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100  m  are  usually  not  important  since  sound  penetrating  more  deeply  will 
be  very  strongly  attenuated.  Routine  coring  to  these  depths  cannot  be 
done,  so  acoustical  (or  seismological )  probes  are  an  essential  part  of 
sediment  measurements. 

There  are  two  general  deep  water  acoustic  propagation  issues  for 

this  program.  First,  the  physical  mechanism  of  sound  attenuation  in 

sediments  needs  to  be.  established  by  laboratory  measurements  and 

analytical  work.  This  effort  will  strengthen  present  capabilities  to 

interpret  and  extrapolate  acoustical  measurements  made  at  sea,  in  both 

region  and  frequency.  A  second  and  distinct  issue  is  the  presence  of 

gas  hydrates  in  marine  sediments.  The  survey  paper  by  Daniels  and 
2 

Vidmar  reviews  what  is  known  about  these  formations,  which  apparently 
are  common  in  the  ocean,  with  particular  attention  to  their  acoustical 
implications. 

In  shallow  water,  acoustic  bottom  interaction  is  a  more  dominant 
factor  than  it  is  in  deep  water.  That  is  due  to  the  fact  that  except  in 
isothermal  water  such  as  in  the  Arctic,  shallow  water  sound  repeatedly 

interacts  with  the  bottom.  However,  the  depth  scales  of  interest  have 
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changed  relative  to  deep  water.  The  work  of  Mitchell  and  Focke,  for 
example,  shows  that  in  shallow  water  the  bottom  properties  are  important 
at  50  Hz  to  50  m  of  sediment,  while  only  a  few  meters  are  important  at 
several  hundred  hertz.  Thus,  the  sediment  properties  important  to  low 
frequency  acoustic  propagation  may  potentially  be  directly  determined  by 
coring. 

As  with  the  deep  water  case,  an  important  shallow  water  acoustic 
propagation  issue  is  the  interpretation  and  extrapolation  of  bottom 
acoustical  measurements.  In  addition,  shallow  water  is  more  likely  to 
contain  coarsely  grained  sand  for  which  some  predictions  of  the  Biot- 
Stoll  theory  can  be  consequential.  Thus,  extension  of  that  theory  to 
real  (variable  grain  size,  shape,  and  composition)  is  important. 


B.  Comments  on  the  Status  of  Theoretical  Modeling  of  Acoustic  Waves 

in  Saturated  Sediments 

An  important  part  of  the  analysis  of  acoustic  waves  in  the  ocean  is 
the  prediction  of  wave  interactions  with  the  bottom.  This  is  a  most 
complex  problem,  not  only  because  of  the  irregular  topography  and  Lhe 
inhomogeneity  of  the  ocean  bottom,  but  because  the  material  comprising 
the  bottom  is  a  fluid  saturated  porous  material  having  very  complicated 
mechanical  characteristics. 

In  this  section,  a  summary  is  given  of  the  theoretical  models  which 
have  been  developed  for  the  analysis  of  wave  propagation  in  saturated 
sediments.  First,  theories  for  wave  propagation  in  the  unbounded  medium 
are  considered.  The  single  component  approach  is  treated  briefly, 
followed  by  a  more  extensive  discussion  of  the  mixture  theory,  due  to 
Biot.  Next,  some  recent  progress  on  the  reflection  and  refraction  of 
waves  at  a  plane  interface  between  the  water  and  a  water  saturated 
sediment  is  described. 

The  simplest  approach  to  analyzing  waves  in  saturated  sediments  is 

to  model  the  saturated  sediment  as  a  single  continuous  medium.  The 
mechanical  properties  of  a  continuous  medium  are  chosen  so  that  its 
behavior  is  equivalent,  in  a  necessarily  limited  sense,  to  that  of  the 

saturated  sediment.  Since  waves  in  saturated  sediments  in  general 

exhibit  dispersion  and  attenuation,  the  equivalent  continuous  medium 
must  be  viscoelastic. 

4-8 

Models  of  this  type  have  been  discussed  by  Hamilton.  He 

suggests  the  use  of  the  particularly  simple  linearly  viscoelastic  model 
in  which  the  Lame  constants  X  and  u  of  an  isotropic,  linearly  elastic 
material  are  replaced  by  complex  moduli  x+ix',  4+iu'.  The  moduli  x, 
u,  and  m’  are  assumed  to  be  frequency  dependent  and  are  to  be  determined 

Q 

from  experiment.  This  type  of  model  is  equivalent  to  a  Voigt  material. 
Hamilton  cites  developments  of  the  model  by  Ferry*®  and  White.** 
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Such  a  purely  empirical  approach  has  substantial  advantages  and 
disadvantages.  The  primary  advantage  is  its  simplicity.  Once  the  coef¬ 
ficients  in  t>e  theory  have  been  evaluated  as  functions  of  frequency, 
then  the  large  body  of  literature  and  many  existing  computer  programs  on 
wave  propagation  in  viscoelastic  materials  can  be  used  to  obtain 
solutions  to  specific  problems.  On  the  other  hand,  since  the  model  does 
not  explicitly  incorporate  the  physical  properties  of  the  saturated 
sediment,  the  coefficients  in  the  theory  must  be  measured  anew  for  each 
change  in  the  properties  of  the  material.  More  importantly,  the  theory 
does  not  have  the  ability  to  predict  the  effects  on  wave  propagation  of 
changes  in  the  material  properties. 

In  order  to  introduce  the  physical  properties  of  the  saturated 
sediment,  the  material  can  be  modeled  as  a  binary  mixture  of  fluid  and 
solid  constituents.  By  retaining  the  individual  identities  of  the  fluid 
and  solid,  it  is  possible  to  express  at  least  part  of  the  coefficients 
in  terms  of  the  physical  properties  of  the  two  constituents. 

A  theory  of  this  type  was  developed  for  application  to  fluid 

12-15 

saturated,  porous  elastic  media  by  Biot.  While  he  derived  his 

equations  on  a  somewhat  intuitive  basis,  many  subsequent  studies  have 
confirmed  that  they  are  consistent  with  the  principles  of  continuum 
mechanics  (see,  for  example,  Bowen, ^  and  Bedford  and  Drumhel ler^) . 

Biot's  equations  can  be  written 
<1-*>  "s  u'(s)m  = 

*  A  e(s)kk (m  *  2“b  e(s)mk,k  *  Q  e(f)kk,ir,  '  b[Cl(s)m-|1(f)m]  '  (1 

*”f  V)m  ■  m  (y ( s)m-u"(  f  )m] 

*  5  e(s)kk,m  *  R  e(v)kk,m  +  b[°( s )m‘Cl{ f )m]  ’  (2' 
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where  <i>  is  the  porosity,  p$  and  are  the  mass  densities  of  the  solid 
and  fluid,  u^m  and  u^m  are  the  displacement  vectors  of  the  solid  and 
fluid,  the  notation  ,m  denotes  partial  differentiation  with  respect  to 
the  coordinate  xm,  e^mk,  and  e^mk  are  the  linear  strains  of  the 
solid  and  fluid  e{s)mk=l/2[u{s)mjk+u(s)kjtT1]  ,  and  A,  «b,  Q,  r,  c,  and  b 
are  constitutive  coefficients. 


Equations  (1)  and  (2)  are  the  equations  of  motion  for  the  solid  and 

fluid  constituents,  respectively.  The  first  term  on  the  right  of  each 

equation  is  a  virtual  mass  term  which  is  linear  in  the  relative 

acceleration  of  the  constituents.  The  last  term  on  the  right  is  a  drag 

term  which  is  linear  in  the  relative  velocity  of  the  constituents.  By 

using  the  solution  for  an  oscillating  cylinder  containing  a  viscous 

fluid,  Biot  was  able  to  evaluate  the  coefficient  of  the  drag  term  b  as  a 

1 3 

function  of  frequency. 

1  O 

In  a  further  important  development,  Biot  and  Willis  showed  that 
the  coefficients  A,  Q,  and  R  could  be  expressed  in  terms  of  the  bulk 
modulus  of  the  solid  material,  the  bulk  modulus  of  the  fluid,  and  the 
bulk  and  shear  moduli  of  the  drained  porous  solid  and 

The  Biot  theory  was  first  applied  to  marine  sediments  by  Stoll  and 
19  ?0-°2 

Bryan  and  Stoll.  u  In  order  to  account  for  dissipative  effects 

associated  with  motion  of  the  granular  matrix  material,  these  authors 

assumed  that  the  moduli  «b  and  ub  were  viscoelastic.  Based  upon 

experimental  data  on  dry  granular  media,  it  was  assumed  that  Kb  and 

were  complex  constants.  It  was  shown  by  Stoll  that  the  theory  could 

correctly  predict  the  variation  with  frequency  of  the  attenuation  of 

on 

compressional  waves  in  saturated  sands. 

The  work  of  Stoll  and  Bryan  was  extended  by  Hovem  and  Ingram,  who 
showed  that  Biot's  technique  for  evaluating  the  drag  coefficient  b  could 
be  used  to  evaluate  the  virtual  mass  coefficient  c  as  well.  They  were 
also  able  to  obtain  a  more  explicit  expression  for  b  in  the  case  of  a 
sediment  consisting  of  particles  of  uniform  size.  They  then  presented 
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very  favorable  comparisons  of  the  theory  with  measurements  of 
congressional  wave  attenuation  in  saturated  sands.  In  addition,  they 
presented  new  experimental  data  on  compressional  waves  in  a  model 
sediment  consisting  of  spherical  glass  beads  saturated  by  water  and 
showed  that  the  theory  correctly  predicted  both  the  attenuation  and 
phase  velocity  of  compressional  waves. 

Thus,  the  extended  Biot  theory  has  been  shown  to  give  quite 
accurate  predictions  in  comparisons  with  direct  measurements  of  both 
phase  velocity  and  attenuation  in  saturated  sediments. 

20 

Stoll  has  used  the'  Biot  theory  to  show  that  the  dominant  loss 
mechanism  at  high  frequencies  is  the  drag  between  the  solid  and  fluid 
constituents,  while  the  dominant  mechanism  at  low  frequencies  is 
dissipation  in  the  viscoelastic  granular  matrix. 

Because  the  scale  of  measurements  of  phase  velocity  and  attenuation 

is  necessarily  small  in  a  laboratory  setting,  the  data  are  in  the 

kilohertz  frequency  range  and  above.  Such  measurements  cannot  be  used 

?4 

to  verify  the  theory  at  lower  frequencies.  Stoll  has  obtained  some 
low  frequency  data  on  the  attenuation  of  shear  waves  in  saturated 
sediments  by  using  a  resonant  column  technique,  but  he  gave  only 
qualitative  comparisons  with  the  theory. 

There  is  a  pressing  need  for  additional  data  at  low  frequencies  and 
for  data  which  can  shed  light  on  the  viscoelastic  properties  of  the 
granular  matrix. 

25 

Recently,  Shirley  et  al.  presented  direct  measurements  of  phase 
velocity  and  attenuation  in  a  glass  bead  sediment  saturated  by  a  mixture 
of  water  and  glycerol.  By  varying  the  proportion  of  glycerol,  the 
viscosity  of  the  liquid  could  be  changed.  In  comparison  with  the 
extended  Biot  theory,  the  data  predicted  a  substantially  greater 
increase  in  attenuation  with  increasing  viscosity.  The  measurements  were 
made  at  a  single  frequency.  Measurements  made  for  a  range  of 
frequencies  would  greatly  aid  in  characterizing  the  viscoelastic 


response  of  the  granular  matrix.  This,  in  turn,  could  conceivably 

permit  extrapolation  of  the  theory  to  low  frequency. 

The  ultimate  objective  of  efforts  to  model  wave  propagation  in 
saturated  sediments  is  the  analysis  of  wave  interactions  with  the  ocean 
bottom.  Although  much  work  remains  to  be  done  on  propagation  in  the 
unbounded  medium,  important  progress  on  the  interaction  problem  has 
recently  been  made  by  Stoll  and  Kan.  They  considered  a  liquid  half¬ 
space  above  a  half-space  of  saturated  sediment,  and  solved  the  problem 
of  the  reflection  and  refraction  of  the  plane  waves  at  the  interface 

using  the  extended  Biot  theory.  The  analysis  was  based  on  the  earlier 

27 

work  on  saturated  porous  media  by  Deresiewicz  and  Rice. 

In  contrast  to  the  results  for  elastic  media,  they  showed  that  the 
reflection  coefficient  exhibited  strong  frequency  dependence,  and  was 
also  a  function  of  the  assumed  value  of  permeability  of  the  sediment. 
The  results  obviously  have  important  implications  for  the  modeling  of 
bottom  interactions. 

The  analysis  of  Stoll  and  Kan  will  have  to  be  extended  before 

realistic  estimates  of  bottom  loss  can  be  achieved.  They  assumed  that 

the  sediment  properties  were  homogeneous  and  isotropic.  It  is,  of 

course,  well  known  that  sediment  properties  are  strongly  depth 
8  22  _ 

dependent.  •  There  is  also  some  evidence  that  sediments  are 
28  29 

anisotropic.  *  Both  of  these  factors  would  influence  the  reflection 
and  refraction  of  waves. 


III.  NEW  LABORATORY  MEASUREMENTS  OF  WAVE  PROPAGATION  DEPENDENCE 

UPON  PORE  FLUID  PROPERTIES 


One  set  of  measurements  of  the  acoustic  properties  of  a  glass  bead 

on 

sediment  with  varied  fluid  makeup  was  taken  in  1980.  This  work  was 
analyzed  and  prepared  for  publication  in  1981;  a  preprint  is  contained 
in  Appendix  A  of  this  report.  The  measurements  were  of  compressional 
wave  and  shear  wave  velocities  and  attenuations  in  a  sediment  with  a 
pore  fluid.  The  content  of  the  fluid  was  varied  from  pure  water  to  28% 
glycerol  in  water.  (The  terms  glycerol  and  glycerine  are  both  used  to 
designate  the  compound  CgHgO^;  the  former  is  preferred  in  the  chemistry 
literature.)  It  was  found  that  the  velocity  measurements  generally 
conformed  to  the  predictions  of  the  Biot-Stoll  theory.  However,  the 
attenuation  data  were  more  scattered,  and  approximately  doubled  with  a 
change  from  0%  to  28%  glycerol.  This  was  in  contrast  to  calculated  the¬ 
oretical  changes  of  attenuations,  which  amounted  to  only  a  few  percent. 
Thus,  it  was  decided  to  repeat  the  measurements  and  try  to  minimize 
sediment  disturbance.  Signal  levels,  the  measured  quantities  which  are 
converted  to  attenuation,  are  most  sensitive  to  sediment  disturbance. 

A.  Physical  Properties 

Variations  in  the  viscosity,  density,  and  bulk  modulus  of  the 
aqueous  solution  were  determined  experimentally,  and  are  shown  in 
Figs.  1  and  2  as  a  function  of  the  percent  concentration  of  glycerol  by 
weight  in  the  solution.  The  viscosity  and  density  properties  were 
checked  against  published  data  for  accuracy. 

The  viscosity  n  (Poise)  was  measured  at  23°  with  an  Ostwalt 

3 

viscometer,  and  the  density  (g/cm  )  with  a  calibrated  100  ml 
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FIGURE  1 

DENSITY  AND  VISCOSITY  versus  PERCENT  WEIGHT  OF  GLYCEROL  IN  SOLUTION 


FIGURE  2 

BULK  MODULUS  versus  PERCENT  WEIGHT  OF  GLYCEROL  IN  SOLUTION 
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picnometer.  The  bulk  modulus  (dyn/cm  )  was  determined  by  measuring 
the  compressional  wave  velocity  and  attenuation  in  a  solution  containing 
various  concentrations  of  glycerol,  then  using  this  data  in  conjunction 
with  the  data  on  density  to  compute  the  bulk  modulus  of  the  solution. 

The  physical  data  that  constitutes  the  characteristics  of  the  glass 
beads  are  shown  in  Table  I  and  were  supplied  by  the  manufacturer  of  the 
beads.  A  different  size  of  beads  than  was  used  in  the  work  of  Appendix  A 
was  selected. 

B.  Preparation  of  the  Sediment 


The  sediment  was  prepared  in  a  steel  tank  (16  cm  x  30.5  cm  x  20  cm 
deep)  by  adding  deionized  water  to  initially  dry  glass  beads,  occupying 
85%  of  the  tank's  volume,  until  the  water  exceeded  the  top  surface  of 
the  beams  by  1  cm.  During  this  procedure,  air  bubbles  are  always 
entrapped  in  the  beads.  To  remove  this  air,  the  tank  was  subjected  to  a 
vacuum  for  24  hours  and  then  heated  to  boiling,  for  about  one  hour.  The 
tank  was  then  allowed  to  cool  to  room  temperature  (approximately  23°C) 
and  then  stirred  with  a  slender  rod  to  visually  inspect  for  any  air 
bubbles  rising  to  the  water's  surface.  After  all  trapped  air  was 
removed,  the  tank  then  became  one  part  of  a  fluid  circulating  system. 
This  apparatus,  shown  schematically  in  Fig.  3,  facilitated  the  variation 
of  the  glycerol  concentration  without  disturbing  the  sediment.  The 
acoustic  transducers  were  embedded  in  the  sediment  to  a  depth  of 
approximately  8  cm,  and  then  the  apparatus  remained  undisturbed  for  at 
least  24  hours.  The  transducers  that  were  used  for  acoustical  measure¬ 
ments  were  similar  to  those  used  for  the  work  described  in  Appendix  A, 
and  consisted  of  a  shear  wave  bender  element  mounted  so  that  the  plane 
of  the  bender  was  vertical.  A  small  compressional  wave  element  was  near 
the  bender  element.  Compressional  wave  data  were  measured  at  a 
frequency  of  114  kHz  and  shear  wave  data  were  measured  at  2.8  kHz. 


TABLE  I 

PHYSICAL  CHARACTERISTICS  OF  SPHERICAL  GLASS  BEADS  OF  SEDIMENT 

Diameter  0.0707  cm 

Density  2.5  gm/cm 

Bulk  modulus  3.5  X  10^  dyn/cm^ 


Porosity 


0.373 


FIGURE  3 

SCHEMATIC  DIAGRAM  OF  EXPERIMENT 
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C.  Acoustical  Measurements 


The  procedure  for  extracting  acoustical  data  at  any  particular 
concentration  of  glycerol  in  the  saturated  fluid  was  to  record  data, 
then  allow  the  apparatus  to  sit  undisturbed  for  24  hours,  after  which 
the  measurements  were  repeated  to  ensure  that  a  stable  sediment  had  been 
obtained.  Once  a  set  of  data  was  taken,  a  new  concentration  of  glycerol 
was  obtained  by  adding  a  quantity  of  pure  glycerol  (usually  400  ml)  to 
the  mixture  tank.  The  fluid  was  then  allowed  to  circulate  in  a  closed 
fashion  by  flowing  through  the  sediment  tank  to  the  reservoir,  from 
which  it  was  pumped  up  to  the  mixture  tank,  creating  a  continual  flow. 
The  viscosity  and  density  of  the  fluid  "ere  tested  periodically  while 
circulation  was  in  progress.  When  successive  tests  indicated  that  these 
properties  were  not  changing  as  circulation  continued,  it  was  decided 
that  a  new  stable  solution  had  been  reached.  The  last  viscosity  and 
density  measurements  recorded  for  this  new  solution  were  then  used  to 
determine  the  concentration  of  glycerol  by  weight  in  the  solution. 
Afterwards,  the  apparatus  remained  undisturbed  for  at  least  24  hours 
before  acoustical  data  for  the  new  solution  were  recorded.  Then  the 
cycle  for  a  new  solution  began. 

Successive  measurements  were  made  to  a  concentration  of  glycerol  by 
weight  of  approximately  35%.  Figures  4-7  give  the  measured  acoustical 
data  plotted  as  a  function  of  glycerol  concentration.  Because  the  bead 
sizes  differed  from  those  of  the  work  reported  in  Appendix  A,  the 
acoustical  properties  are  different.  This  is  clearly  evident  when  shear 
wave  velocities  are  compared. 

Theoretical  values  of  wave  velocities  and  attenuations  have  been 
computed  and  are  plotted  as  solid  curves  on  Figs.  4-7.  Although  the 
diameter  of  the  glass  beads  used  in  obtaining  the  present  data  was 
approximately  four  times  the  diameter  of  the  beads  that  were  used  in  our 
earlier  work,  several  significant  qualitative  comparisons  can  be  made. 
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COMPRESSION AL  WAVI 


FIGURE  4 

COMPRESSIONAL  WAVE  VELOCITY  IN  THE  SATURATED  SEDIMENT 
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FIGURE  5 

COMPRESSION AL  WAVE  ATTENUATION  IN  THE  SATURATED  SEDIMENT 
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FIGURE  6 

SHEAR  WAVE  VELOCITY  IN  THE  SATURATED  SEDIMENT 
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FIGURE  7 

SHEAR  WAVE  ATTENUATION  IN  THE  SATURATED  SEDIMENT 
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The  agreement  between  the  measured  and  theoretical  congressional 
wave  velocities  (Fig.  4)  is  approximately  the  same  as  that  of  the 
earlier  work.  In  the  earlier  work,  the  data  on  the  compressional  and 
shear  wave  attenuations  showed  larger  increases  with  increasing  glycerol 
concentration  than  were  predicted  by  the  theory.  This  is  not  observed 
in  the  new  results  (Figs.  5  and  7).  Also,  in  the  earlier  work  the  data 
on  the  shear  wave  velocity  were  found  to  decrease  more  rapidly  than  was 
predicted  by  the  theory.  In  the  new  results,  the  data  actually 
increase,  then  decrease  with  increasing  glycerol  concentration — an 
observation  which  cannot  be  accounted  for  by  the  theory. 

Thus,  the  new  data  are  more  consistent  with  the  theoretical 
predictions  of  the  attenuation,  which  may  be  attributable  to  the 
improved  experimental  procedure  that  was  used.  However,  the  new  data 
also  exhibit  anomalous  features  (Fig.  6)  which  require  further 
investigation. 
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Wave  propagation  in  a  saturated  model  sediment  with  varied  liquid  properties 

D.  J.  Shirley3^,  A.  Bedford^  and  S.  K.  Mitchell 

Applied  Research  Laboratories,  The  University  of  Texas  at  Austin, 
Austin,  Texas  78712 

Experimental  measurements  have  been  made  of  the  phase  velocity  and 
attenuation  of  compressional  and  shear  waves  in  a  bed  of  glass  beads 
saturated  by  a  mixture  of  water  and  glycerine.  By  varying  the 
proportion  of  glycerine,  a  ranqe  of  values  of  the  viscosity,  density 
and  bulk  modulus  of  the  mixture  was  obtained.  The  data  have  been 
compared  with  theoretical  predictions  based  upon  the  Biot  theory. 

INTRODUCTION 

In  order  to  predict  the  effects  of  the  interaction  of  acoustic  waves 

with  the  ocean  bottom,  it  is  necessary  to  have  a  comprehensive  understanding 

of  the  propagation  of  waves  in  saturated  sediments.  Great  progress  toward 

1  2 

such  an  understanding  has  recently  been  achieved  by  Stoll  and  Bryan  ,  Stoll 

5 

and  Hovem  and  Ingram  .  These  authors  analyzed  the  propagation  of  waves  in 

6-9 

sediments  by  using  the  theory  of  fluid  saturated  porous  media  due  to  Biot 
Biot  modeled  the  porous  medium  and  the  fluid  as  superimposed  continua, 
and  used  a  variational  approach  to  derive  equations  of  motion  for  the 
individual  constituents.  The  equations  contained  drag  terms  associated  with 
the  relative  velocity  of  the  constituents,  and  also  virtual  mass  terms 
associated  with  their  relative  acceleration.  By  using  the  solution  for  a 
cylinder  containing  viscous  fluid  and  subjected  to  axial  oscillations,  Biot 


was  able  to  evaluate  the  coefficient  of  the  drag  term  as  a  function  of 

frequency.  Furthermore,  Biot  and  Willis^  showed  that  the  elastic  coefficients 

in  Biot's  equations  could  be  expressed  in  terms  of  the  bulk  moduli  of  the 

fluid  and  solid  materials  and  the  elastic  moduli  of  the  porous  solid.  In 

this  regard  also  see  the  discussion  by  Stoll^. 

Stoll  and  Bryan  were  the  first  to  apply  Biot's  equations  to  marine 

sediments.  Unlike  the  porous  elastic  material  which  Biot  had  considered, 

the  granular  matrix  materials  in  sediments  were  known  to  exhibit  dissipative 

effects.  They  therefore  assumed  that  the  moduli  of  the  solid  material  were 

complex.  On  the  basis  of  experimental  data  on  dry  granular  materials,  they 

1 2 

assumed  that  these  complex  moduli  were  not  frequency  dependent  .  Using 

this  approach,  Stoll  subsequently  showed  that  the  theory  could  correctly 

predict  the  variation  with  frequency  of  the  attenuation  of  compressional 

1 3 

waves  in  saturated  sands  .  He  also  pointed  out  that  dissipation  due  to 
relative  motion  of  the  fluid  and  solid  was  the  dominant  loss  mechanism  at 
higher  frequencies,  while  losses  associated  with  the  granular  matrix  were 
dominant  at  low  frequency. 

Hovem  and  Ingram  were  able  to  show  that  the  technique  used  by  Biot  to 

evaluate  the  drag  term  in  the  equations  could  also  be  used  to  evaluate  the 

14 

coefficient  of  the  virtual  mass  term  as  a  function  of  frequency  .  They 
showed  that  the  theory  compared  favorably  with  measurements  of  compressional 
wave  attenuation  in  saturated  sands.  In  addition,  they  presented  new 
experimental  data  on  compressional  waves  in  a  model  sediment  consisting  of 
spherical  glass  beads  saturated  by  water,  and  showed  that  the  theory  correctly 
predicted  both  the  attenuation  and  phase  velocity  as  functions  of  frequency. 


In  this  paper,  results  are  presented  of  measurements  of  phase  velocity 
and  attenuation  of  congressional  and  shear  waves  in  a  glass  bead  model 
sediment  of  the  type  used  by  Hovem  and  Ingram.  The  novel  aspect  of  the 
measurements  was  that  a  mixture  of  water  and  glycerine  was  used  as  the 
saturating  fluid.  By  changing  the  proportion  of  glycerine  in  the  fluid, 
data  was  obtained  for  a  range  of  values  of  the  viscosity,  density  and  bulk 
modulus  of  the  fluid.  Since  these  parameters  appear  explicitly  in  the 
extended  Biot  equations  developed  by  Stoll -Bryan  and  Hovem-Ingram,  the 
measurements  provide  a  critical  test  of  the  theory. 

I.  EXPERIMENTAL  MEASUREMENTS 

When  glycerine  is  added  to  water,  the  density,  viscosity  and  bulk  modulus 

of  the  mixture  are  larger  than  ‘.he  values  for  water.  Published  data  are 

available  for  the  density  and  viscosity  of  a  water-qlycerine  solution  at 

1 5 

20°C  as  a  function  of  the  percent  of  glycerine  by  weight  .  For  0  £  W  £  28  , 
where  W  is  the  percent  e*  glycerine  by  weight,  the  data  on  the  density 
Pf  (g/cm  )  and  viscosit/  n  (Poise)  of  the  mixture  can  be  approximated  by 
the  empirical  equation' 

pf  =  0.99-12  +  (2.284  x  10"3)  W 

+  (4.067  x  10'6)  W2  +  (4.960  x  10"8)  W3  ,  (1) 

n  -  0.01002  +  (2.461  x  10"4)  W 


As  part  of  this  study,  measurements  were  made  of  the  wave  velocity  and 
attenuation  in  the  water-glycerine  mixture.  The  measurements  were  made  at  a 
frequency  of  114  kHz  (the  higher  of  the  two  frequencies  used  in  the  model 
sediment  measurements).  The  wave  velocity  data  are  shown  in  Fig.  1.  No 
variation  in  the  wave  attenuation  due  to  the  addition  of  glycerine  was 
observed.  Once  the  wave  velocity  was  known,  it  was  used  together  with  the 
data  on  the  density  to  compute  the  bulk  modulus  of  the  mixture  as  a  function 
of  the  percent  of  glycerine  by  weight  (Fig.  2).  For  0  ^  W  ^  28  ,  the  bulk 

p 

modulus  of  the  mixture,  (dvne/cm  )  ,  can  be  approximated  by 

Kf  =  [2.175  +  (2.052  x  10"2)  W 

+  (1.400  x  10"4)  W2  -  (1.667  x  10-7)  W3]  x  1010  .  (3) 

The  solid  matrix  material  consisted  of  a  bed  of  nominally  spherical  glass 

O 

beads  of  0.0177  cm  diameter.  The  density  of  the  glass  was  2.5  g/cm  and  the 
porosity  of  the  bed  was  0.365. 

The  measurements  of  wave  velocity  and  attenuation  were  made  in  a 
sediment  tank  16  cm  x  30  cm  x  20  cm  deep.  The  depth  of  the  transducers  in 
the  sediment  was  approximately  10  cm. 

The  sediment  sample  was  prepared  by  adding  demineralized  water  to  the 
initially  dry  material,  boiling  the  mixture  and  then  subjecting  the  cooled 
sediment  to  a  vacuum  for  24  hours.  Once  the  sediment  was  ready  for  measurement, 
the  transducers  were  inserted  into  the  material  and  the  apparatus  allowed  to 
remain  undisturbed  for  another  24  hours.  Acoustical  measurements  were  then 
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made  and  the  sediment  was  allowed  to  sit  undisturbed  for  another  24  hours, 
after  which  the  acoustical  measurements  were  repeated.  This  procedure  was 
repeated  until  successive  shear  wave  measurements  we re  essentially  identical. 

It  was  usually  found  that  the  sediment  was  stabilized  by  the  third  or 
fourth  24  hour  interval. 

Once  a  set  of  measurements  had  been  made,  the  sediment  was  removed  from 
the  tank,  approximately  400  ml  of  glycerin  added  and  the  new  mixture  thoroughly 
stirred.  The  sediment  was  again  evacuated  to  remove  entrained  air  and  the 
measurement  procedure  described  above  was  repeated.  After  the  acoustical 
measurements  were  completed  for  a  particular  concentration  of  glycerin,  a 
sample  of  the  pore  fluid  was  removed  from  the  tank  for  viscosity  and  density 
measurements.  Viscosity  was  measured  at  20°C  with  a  modified  Ostwalt 
viscometer  and  density  was  measured  at  20°C  with  a  calibrated  50  ml  picnometer. 
The  viscosity  and  density  measurements  were  used  to  determine  concentration. 

After  successive  measurements  to  a  concentration  of  approximately  25% 
glycerin,  the  sediment  was  discarded  and  the  whole  procedure  repeated  with 
freshly  prepared  sediment  as  a  check. 

The  transducers  that  were  used  to  make  the  compressional  and  shear  wave 
measurements  were  similar  to  those  described  by  Shirley16,  and  consisted 
of  a  shear  wave  bender  element  mounted  so  that  the  plane  of  the  bender  was 
vertical,  and  a  small  compressional  wave  element  near  the  bender  element. 

One  projector  was  used  with  two  receivers  so  that  attenuation  could  be 
calculated  from  the  difference  in  amplitude  of  signals  from  the  receivers. 
Compressional  v/ave  data  were  measured  at  a  frequency  of  114  kHz  and  shear 
wave  data  were  measured  at  2.8  kHz. 
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In  Fig.  3,  the  measured  compressional  phase  velocity  is  plotted  as  a 
function  of  glycerine  concentration.  (The  curves  on  Figs.  3-6  will  be 
discussed  in  Part  II  of  the  paper.)  The  phase  velocity  increases  approximately 
linearly  with  increasing  concentration  and  there  is  little  scatter  in  the 
data.  The  increase  in  the  velocity  occurs  primarily  due  to  the  increasing 
bulk  modulus  of  the  fluid  (Fig.  2). 

Note  that  the  measured  value  of  the  phase  velocity  with  no  glycerine 
present  differs  from  the  compressional  phase  velocity  measured  at  the  same 
frequency  by  Hovem  and  Ingram^  by  over  100  m/s.  The  difference  is  due  to 
Hovem  and  Ingram's  measurements  being  made  at  a  depth  of  20  cm  in  the 
sediment,  compared  to  a  depth  of  10  cm  for  the  present  measurements.  This 

is  an  example  of  the  strong  depth  dependence  of  the  acoustic  properties  of 

* 

18  19 

saturated  sediments  *  . 

The  measured  compressional  wave  attenuation  is  plotted  in  Fig.  4.  The 
data  exhibits  substantial  scatter,  but  suggests  an  approximately  linear 
increase  in  attenuation  with  increasing  concentration.  Increased  attenuation 
would  be  expected  since  the  drag  arising  from  the  relative  motion  of  the 
fluid  and  solid  increases  with  increasing  viscosity  of  the  fluid.  However, 
it  will  be  shown  in  Part  II  that  this  mechanism  alone  does  not  explain  the 
observed  increase. 

The  measured  shear  phase  velocity  is  plotted  in  Fig.  5.  Again  there 
is  substantial  scatter,  but  the  data  suggests  a  decrease  in  the  velocity  with 
increasing  concentration. 

The  measured  shear  wave  attenuation  is  plotted  in  Fig.  6.  As  with  the 
compressional  attenuation,  the  data  indicate  a  substantial  increase  in 
attenuation  with  increasing  glycerine  concentration. 
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II.  THEORETICAL  COMPARISONS 

20 

Biot's  equations  can  be  written,  in  Cartesian  tensor  notation, 

32  32 

(1  -  *)  y(s),'  -  c  ^5-  (u(s)m  -  u(f)m> 

+  Ae(s)kk,m  +  ^ybe(s)mk,k 

+  ^  e(f)kk,m 

-  b  It  (u(s)m  '  u(f)n.>  •  <4> 

2  2 

^  Pf  at2  U(f)m  =  C  at2  (U(s)m  "  U(f)m) 

+  ^  e(s)kk,m  +  R  e(f)kk,m 

+  b  It  ^u(s)m  "  U(f)m^  *  ^ 

where  4>  is  the  porosity,  p$  and  are  the  mass  densities  of  the  solid 
and  fluid,  u(s)m  and  u(f)m  are  tbe  displacement  vectors  of  the  solid  and 
fluid,  the  notation  ,m  denotes  partial  differentiation  with  respect  to  the 
coordinate  xm  ,  e(s)m|<  and  are  the  linear  strains  of  the  solid 

and  fluid  (e(s)m|(  =■  'j  (u(s)m>k  *  u(s)k>m))  ,  and  A  ,  Jb  ,  Q  ,  R  ,  c  and 
b  are  constitutive  coefficients.  The  first  terms  on  the  right  of  Eqs.  (4) 
and  (5)  are  the  virtual  mass  terms  and  the  last  terms  on  the  right  are  the 
drag  terms. 
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By  introducing  the  notation 


u(s)m,m  5  ^  ^  ^u(s)m,m  "  u(f)m,n^ 


and  taking  the  divergence  of  Eqs.  (4)  and  (5),  they  can  be  written 


T  <P<=  -  ofc)  -  He>mm  -  Ccl|njn  , 


(ofe  -  pc0  -  Ce,mm  -  Mc,r 


+  b  l_  r 

+  "7  3t  - 

<P 


where  p  =  (1  -  <J>)  p$  +  $  p^.  is  the  mass  density  of  the  fluid-solid  mixture 


H  =  A  +  2yb  +  2Q  +  R  , 


C  =  (Q  +  R)  /  <t>  , 


M  =  R  /  <J> 


Pc  s  Pf  /  <p  +  c  /  <p  . 
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Equations  (7)  and  (8)  are  essentially  the  forms  of  Biot's  equations  for 

congressional  waves  used  by  Stoll  and  Hovem  and  Ingram.  Hovem  and  Ingram 

21 

extended  Biot's  analysis  and  evaluated  the  coefficients  b  and  c  as 
functions  of  frequency.  Their  results  were 


b  /  <t>  =  n  Real  {F{k))  /  BQ  , 


c  /  <j>  =  n  Imaginary  (F(<))  /  (coBQ)  , 


where  n  is  the  viscosity  of  the  liquid,  u  is  the  frequency  (rad/s),  and 


<  ~  ap  (topf  /  n)‘ 


The  term  ap  is  a  pore  size  parameter  defined  by 


ap  =  dm  I  [3(1  -  *)] 


and  dm  is  the  diameter  of  the  spherical  grains  of  the  sediment.  The  function 
F(k)  is22 


F(k)  =  l  kT(k)  /  [1  -  2T(<)  /  i<]  , 


where  T(k)  is  the  Kelvin  function 


T(k)  =  [ber'(ic)  +  i  bei'(ic)]  /  [ber(K)  +  i  bei (vc) ]  . 

35 


(14) 


jermeability  BQ  is 

Bo  ’  ♦  ap  /  4  *  • 

( 

where  k  is  the  Kozeny-Carman  constant.  For  the  glass  beads  used  in  the 
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present  study.  Bell  obtained  a  value  k  =  3.98  on  the  basis  of  direct 

measurements  of  the  permeability. 

The  expressions  for  the  coefficients  H  ,  C  and  M  in  terms  of  the 
bulk  moduli  of  the  solid  and  fluid  materials  K$  ,  Kf  and  the  bulk  and 
shear  moduli  of  the  drained  granular  matrix  material  i<b  ,  ub  are2^ 

H  =  (Ks  -  i<b)2  /  (0  -  Kb)  +  Kb  +  4ub  /  3  , 

C  =  Ks  (Ks  '  V  '  (D  ’  V  ’ 

( 

M  =  K2  /  (0  -  Kb)  , 

where 

D  »  Ks  [1  +  <J>  (Ks  /  Kf  -  1)]  . 

( 

By  introducing  the  notation 

<j  =  curl  u,  \  ,0  =  4)  curl  (u,_>  -  u,f\) 

( 
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described  in  Part  I.  With  the  exception  of  the  moduli  Kb  and  ub  , 
parameters  in  the  equations  are  known  in  terms  of  the  physical  properties  of 
the  fluid  and  solid  constituents,  and  are  summarized  in  Table  I. 

-  _  2g 

In  determining  and  y^  it  was  assumed,  following  Hovem  and  Ingram  , 

that  =  ub  ,  which  corresponds  to  a  Poisson's  ratio  of  the  solid  matrix 

27 

material  of  0.125.  This  is  within  the  range  of  0.1  to  0.17  which  Stoll 

suggests  for  a  granular  material.  The  values  of  y^  and  6s  were  then 

chosen  so  that  the  theoretical  values  of  the  shear  phase  velocity  and 

attenuation  agreed  with  the  data  for  zero  percent  glycerine.  (Note  from  Eqs. 

(19)  and  (20)  that  the  shear  phase  velocity  and  attenuation  are  independent 

of  «b  .)  It  was  found  that  values  of  ub  =  8  x  10  dyne/cm  and  6s  =  0 

gave  a  shear  phase  velocity  of  66.94  m/s  and  a  shear  attenuation  of  68.41  dB/m. 

Increasing  6$  to  0.1  gave  a  shear  phase  velocity  of  67.00  m/s  and  a  shear 

attenuation  of  104.70  dB/m,  which  agree  with  the  averages  of  the  measured 

values.  Next,  using  the  latter  values  for  y^  and  and  assuming  that 

Kb  =  ub  anc*  6C  =  0  9ave  a  comPress'’ona^  Phase  velocity  of  1807.89  m/s 

and  a  compressional  attenuation  of  30.91  dB/m.  These  values  agree  quite  well 

with  the  measurements.  Increasing  6c  to  0.1  had  a  very  small  effect  on 

the  compressional  attenuation,  increasing  it  to  30.93  dB/m. 

7  2 

Then,  using  the  values  =  y^  =  8  x  10  dyne/cm  and  6s  =  &c  =  0.1  , 
the  -.ompressional  and  shear  phase  velocity  and  attenuation  were  determined 
as  functions  of  glycerine  concentration.  The  results  are  indicated  by  the 
solid  curves  in  Figs.  3-6. 

The  increase  in  the  compressional  phase  velocity  with  increasing 
glycerine  concentration,  due  primarily  to  the  increasing  bulk  modulus  of  the 
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fluid,  is  predicted  very  well  by  the  theory  (Fig.  3).  The  theory  predicts 
a  slight  decrease  in  the  shear  phase  velocity  with  increasing  concentration 
due  to  the  increasing  density  of  the  fluid.  The  data  suggest  a  larger 
decrease,  but  this  remains  conjectural  due  to  the  scatter  of  the  data. 

For  both  compress ional  and  shear  attenuation,  the  data  indicate  a 
substantially  larger  increase  with  increasing  concentration  than  is  predicted 
by  the  theory.  The  slight  increase  in  attenuation  predicted  by  the  theory 
is  due  to  the  dependence  of  the  drag  on  the  viscosity  of  the  fluid.  This 
mechanism  clearly  does  not  account  for  the  measured  attenuation.  Thus  the 
data  indicate  that  there  is  viscous  dissipation  present  which  is  a  function 
of  the  viscosity  of  the  liquid. 

The  possibility  of  viscous  dissipation  of  this  type  has  been  discussed  by 

OO 

Stoll  .  He  stated  that  "...  in  a  general  model,  it  may  be  necessary  to 
allow  for  both  frictional  and  viscous  losses  in  the  parameters  that  describe 
the  response  of  the  skeletal  frame,  particularly  if  high  frequencies  are  to 
be  considered." 

For  the  shear  waves,  a  determination  was  made  of  the  value  of  the  shear 
logarithmic  decrement  necessary  for  the  theory  to  predict  the  measured 
attenuation.  This  was  done  for  each  data  point  in  Fig.  6.  The  resulting 
values  are  plotted  in  Fig.  7  as  a  function  of  the  viscosity  of  the  fluid.  The 
logarithmic  decrement  increases  approximately  linearly  with  increasing 
viscosity  of  the  fluid.  It  is  interesting  to  note  from  Fig.  7  that  the 
data  appear  to  be  approaching  the  horizontal  axis  at  a  finite  value  of  the 
viscosity  (=  0.007  Poise).  There  may  be  a  threshold  value  of  the  viscosity 
below  which  the  fluid  does  not  affect  the  viscoelastic  response  of  the  matrix 
material . 
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The  values  of  the  shear  logarithmic  decrement  plotted  in  Fig.  7  were 
approximated  by  the  linear  relation  (shown  on  the  figure) 

6S  =  -0.241  +  35. 7n  .  (22) 

Equation  (22)  was  then  used  in  the  theory  to  determine  the  shear  wave 
attenuation  as  a  function  of  glycerince  concentration.  The  result  is  shown 
as  the  dashed  curve  in  Fig.  6. 

The  situation  with  regard  to  interpretation  of  the  compressional  wave 
data  was  more  complex.  While  the  shear  wave  attenuation  does  not  depend 
upon  Kb  ,  the  compressional  wave  attenuation  depends  upon  both  Kb  and 
ub  .  Furthermore,  the  compressional  and  shear  measurements  were  made  at 
different  frequencies.  Since  viscous  dissipation  would  depend  upon  frequency, 
the  shear  measurements  could  not  be  used  to  evaluate  pb  at  the  frequency 
at  which  the  compressional  measurements  were  made. 

In  order  to  obtain  an  estimate,  it  was  assumed  that  the  compressional 
and  shear  logarithmic  decrements  were  equal,  6c  =  6$  .  It  was  found  that 
for  26  percent  glycerine,  a  value  6c  =  35  resulted  in  a  compressional 
phase  velocity  of  1983.67  m/s  and  a  compressional  attenuation  of  53.11  dB/m, 
in  agreement  with  the  data.  Note  that  this  logarithmic  decrement  is 
approximately  seventy  times  the  value  of  <5  required  to  match  the  shear 
wave  attenuation.  The  ratio  of  the  frequency  for  the  compressional 
measurements  to  the  frequency  for  the  shear  measurements  was  approximately 
forty.  This  strongly  suggests  that  the  loss  parameters  are  frequency 
dependent. 
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It  is  obvious  that  the  value  35  is  completely  unrealistic  if  6c  is 
interpreted  literally  as  a  logarithmic  decrement.  The  reason  for  the  large 
value  can  be  seen  by  an  examination  of  Eqs.  (7),  (16)  and  (17).  Note  from 
Eq.  (7)  that  H  is  the  compressional  modulus  of  the  fluid-solid  mixture 
when  the  strains  in  the  fluid  and  solid  are  equal.  Since  <<  Ks  in  the 
case  being  considered,  it  can  be  seen  from  Eqs.  (16)  and  (17)  that  a 
relatively  large  change  in  is  necessary  to  cause  a  substantial  change 

I 

in  H  .  The  magnitude  of  should  perhaps  be  compared  to  the  real  part 
of  H  ,  rather  than  with  .  In  the  present  case,  /  Real  (H)  =  0.0132. 
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Figure  Captions 


Fig.  1  Congressional  wave  velocity  in  the  water-glycerine  mixture. 

Fig.  2  Bulk  modulus  of  the  water-glycerine  mixture. 

Fig.  3  Congressional  phase  velocity  in  the  saturated  sediment. 

Fig.  4  Compressional  wave  attenuation  in  the  saturated  sediment. 

Fig.  5  Shear  phase  velocity  in  the  saturated  sediment. 

Fig.  6  Shear  wave  attenuation  in  the  saturated  sediment. 

Fig.  7  Computed  shear  logarithmic  decrement  as  a  function  of  fluid 

viscosity. 
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COMPRESSIONAL  ATTENUATION  -  dB/m 


PERCENT  GLYCERINE  BY  WEIGHT 


SHEAR  ATTENUATION 


Table  I.  Material  Constants  and  Parameters 


Porosity,  <p 
Bead  diameter,  d 

m 

Solid  material  density,  ps 
Fluid  density, 

Fluid  viscosity,  n 
Solid  material  bulk  modulus 
Fluid  bulk  modulus, 
Kozeny-Carman  constant,  k 


0.365 
0.0177  cm 
2.5  g/cm^ 

From  Eq.  (1 ) 

From  Eq.  (2) 

K$  3.5  x  10"  dyne/cm' 

From  Eq.  (3) 

3.98 


50 


^Present  address:  Southwest  Research  Institute,  San  Antonio,  TX  78284. 

^On  leave  from  Department  of  Aerospace  Engineering  and  Engineering  Mechanics, 
The  University  of  Texas  at  Austin,  Austin,  TX  78712. 

^R.  D.  Stoll  and  G.  M.  Bryan,  "Wave  attenuation  in  saturated  sediments," 

J.  Acoust.  Soc.  Am.  47,  1440-1447  (1969). 

^R.  D.  Stoll,  "Acoustic  waves  in  saturated  sediments,"  in  Physics  of  Sound  in 
Marine  Sediments,  edited  by  L.  D.  Hampton  (Plenum,  New  York,  1974JT 

3R.  D.  Stoll,  "Acoustic  waves  in  ocean  sediments,"  Geophys.  42,  715-725  (1977). 

4R.  D.  Stoll,  "Theoretical  aspects  of  sound  transmission  in  sediments,"  J. 
Acoust.  Soc.  Am.  68,  1341-1350  (1980). 

C  * 

30.  M.  Hovem  and  G.  D.  Ingram,  "Viscous  attenuation  of  sound  in  saturated  sand," 
J.  Acoust.  Soc.  Am.  66,  1807-1812  (1979). 

A.  Biot,  "Theory  of  elastic  waves  in  a  fluid  saturated  porous  solid. 

I.  Low  frequency  range,"  J.  Acoust.  Soc.  Am.  28,  168-178  (1956). 

7M.  S.  Biot,  "Theory  of  elastic  waves  in  a  fluid  saturated  porous  solid. 

II.  Higher  frequency  range,"  J.  Acoust.  Soc.  Am.  28,  179-191  (1956). 

8M.  A.  Biot,  "Mechanics  of  deformation  and  acoustic  propagation  in  porous 
media,"  J.  Appl .  Phys.  33,  1482-1498  (1962). 

9M.  A.  Biot,  "Generalized  theory  of  acoustic  propagation  in  porous  dissipative 
media,"  J.  Acoust.  Soc.  Am.  34,  1254-1264  (1962). 

'8M.  A.  Biot  and  D.  G.  Willis,  "The  elastic  coefficients  of  the  theory  of 
consolidation,"  J.  Appl.  Mech.  24,  594-601  (1957). 

nRef.  2,  pp.  21-25. 

12Ref.  1,  1442. 

13Ref.  2,  Fig.  4. 

14Ref.  5,  pp.  1807-1809. 


51 


1C 

CRC  Handbook  of  Chemistry  and  Physics,  edited  by  R.  C.  Weast  (CRC  Press, 

West  Palm  Beach,  horida,  1978),  pp.  D-277,  278. 

^80.  J.  Shirley,  "An  improved  shear  wave  transducer,"  J.  Acoust.  Soc.  Am.  63, 
1643-1645  (1978). 

17Ref.  5,  Fig.  3. 

18Ref.  4,  pp.  1343,  1346. 

1 9 

E.  L.  Hamilton,  "Geoacoustic  model inq  of  the  sea  floor,"  J.  Acoust.  Soc.  Am. 
68,  1313-1340  (1980). 

20Ref.  6,  Eq.  (6.7). 

21Ref.  7,  pp.  180-185. 

22 

It  should  be  noted  that  Ref.  5  contains  misprints  in  the  equations  for 
F(k)  and  T(k). 

23 

D.  W.  Bell,  "Shear  wave  propagation  in  unconsolidated  fluid  saturated  porous 
media,"  Report  ARL-TR-79-31 ,  Applied  Research  Laboratories,  University  of 
Texas  at  Austin  (1979). 

24Ref .  2,  Eq.  (13). 

25Ref.  2,  p.  32. 

26Ref .  5,  p.  1810. 

27Ref.  4,  p.  1345 

28Ref .  4,  pp.  1343-1345 


DOCUMENTATION  UNDER  ARL:UT 
ONR  CODE  425  PROGRAM 


S.  R.  Addy,  E.  W.  Behrens,  T.  R.  Haines,  D.  J.  Shirley,  and  J.  L.  Worzel, 
"Correlation  of  Some  Lithologic  and  Physical  Characteristics  of  Sediments 
with  High  Frequency  Subbottom  Reflection  Types,"  Proceedings  of  the  11th 
Annual  Offshore  Technology  Conference,  Houston,  Texas,  30  April  -  3  May 
1979. 

A.  L.  Anderson,  "Acoustics  of  Gas-Bearing  Sediments,"  Applied  Research 
Laboratories  Technical  Report  No.  74-19  (ARL-TR-74-19) ,  Applied  Research 
Laboratories,  The  University  of  Texas  at  Austin,  May  1974. 

A.  L.  Anderson  and  L.  D.  Hampton,  "In  Situ  Measurement  of  Sediment  Acoustic 
Properties  During  Coring,"  presented  at  the  ONR  Symposium  on  Physical  and 
Engineering  Properties  of  Deep-Sea  Sediments,  Airlie  House,  Airlie,  Vir¬ 
ginia,  24-27  April  1973. 

A.  L.  Anderson  and  L.  D.  Hampton,  "Measurement  of  In  Situ  Acoustic  Proper¬ 
ties  During  Sediment  Coring,"  presented  at  the  ONR  Symposium  on  the  Physics 
of  Sound  in  Marine  Sediments,  Lakeway  Inn,  Austin,  Texas,  8-10  May  1973. 

A.  L.  Anderson  and  L.  D.  Hampton,  “A  Method  for  Measuring  In  Situ  Acoustic 
Properties  During  Sediment  Coring,"  In  Physics  of  Sound  in  Marine  Sediments, 
edited  by  Loyd  Hampton  (Plenum  Press,  1974K 

A.  L.  Anderson  and  L.  D.  Hampton,  "In  Situ  Measurements  of  Sediment  Acous¬ 
tic  Properties  During  Coring,"  In  Deep-Sea  Sediments,  Physical  and  Mechanical 
Properties,  edited  by  A.  L.  Inderbitzen  (Plenum  Press,  1974Y. 

A.  L.  Anderson  and  L.  D.  Hampton,  "Use  of  Tubes  for  Measurement  of  Acous¬ 
tical  Properties  of  Materials,"  presented  at  the  89th  Meeting  of  the  Acous¬ 
tical  Society  of  America,  Austin,  Texas,  April  1975. 

A.  L.  Anderson  and  L.  D.  Hampton,  "Acoustics  of  Gas-Bearing  Sediments. 

Part  I:  Background,"  Applied  Research  Laboratories  Technical  Paper  No. 

79-63  (ARL-TP-79-63),  Applied  Research  Laboratories,  The  University  of 
Texas  at  Austin,  October  1979. 

A.  L.  Anderson  and  L.  D.  Hampton,  "Acoustics  of  Gas-Bearing  Sediments. 

Part  II:  Measurements  and  Models,"  Applied  Research  Laboratories  Technical 
Paper  No.  79-64  (ARL-TP-79-64 ) ,  Applied  Research  Laboratories,  The  Univer¬ 
sity  of  Texas  at  Austin,  October  1979. 

A.  L.  Anderson  and  L.  0.  Hampton,  "Acoustics  of  Gas-Bearing  Sediments. 

Part  I:  Background,"  J.  Acoust.  Soc.  Am.  67_,  1865-1889  (1980). 


A.  L.  Anderson  and  L.  0.  Hampton,  "Acoustics  of  Gas-Bearing  Sediments. 

Part  II:  Measurements  and  Models,"  0.  Acoust.  Soc.  Am.  67^,  1890-1903  (1980). 

A.  L.  Anderson,  R.  J.  Harwood,  and  L.  0.  Hampton,  "Temperature  Studies 
on  Lake  Travis  Stratification  in  a  Warm  Monomictic  Reservoir,"  The  Texas 
Journal  of  Science,  Vol.  XXVI,  Nos.  3  and  4,  353-371,  August  1975. 

D.  W.  Bell,  "Shear  Wave  Propagation  in  Unconsolidated  Fluid  Saturated  Porous 
Media,"  Dissertation  for  Ph.D.  awarded  by  The  University  of  Texas  at  Austin, 
May  1979.  Also  Applied  Research  Laboratories  Technical  Report  No.  79-31 
(ARL-TR-79-31),  Applied  Research  Laboratories,  The  University  of  Texas 
at  Austin,  May  1979. 

D.  W.  Bell  and  0.  J.  Shirley,  "Temperature  Variation  of  the  Acoustical 
Properties  of  Laboratory  Sediments,"  Applied  Research  Laboratories  Tech¬ 
nical  Paper  No.  80-1  (ARL-TP-80-1) ,  Applied  Research  Laboratories,  The 
University  of  Texas  at  Austin,  January  1980. 

D.  W.  Bell  and  D.  J.  Shirley,  "Temperature  Variation  of  the  Acoustical 
Properties  of  Laboratory  Sediments,"  J.  Acoust.  Soc.  Am.  68,  227-231  (1980). 

L.  D.  Hampton,  "ARL  Experience  with  Acoustics  and  Gas  in  Sediments,"  pre¬ 
sented  at  Symposium  on  Natural  Gases  in  Marine  Sediments  and  Their  Mode 
of  Distribution,  The  University  of  California  Lake  Arrowhead  Conference 
Center,  Lake  Arrowhead,  California,  28-30  November  1972. 

L.  D.  Hampton  and  A.  L.  Anderson,  "Acoustics  and  Gas  in  Sediments":  Applied 
Research  Laboratories  (ARL)  Experience,"  in  Natural  Gases  in  Marine  Sediments, 
Marine  Science  Vol.  Ill,  Edited  by  Isaac  R.  Kaplan  (Plenum  Press,  1974). 

J.  M.  Hovem,  "Some  Aspects  of  Sound  Propagation  in  Saturated  Sand,"  presented 
at  the  Seminar  on  Bottom  Effects  in  Underwater  Sound  Propagation,  Miami, 
Florida,  26-28  April  1979. 

J.  M.  Hovem,  "The  Nonlinearity  Parameter  of  Saturated  Marine  Sediments," 
Applied  Research  Laboratories  Technical  Paper  No.  79-17  (ARL-TP-79-17), 

Applied  Research  Laboratories,  The  University  of  Texas  at  Austin,  revised 
May  1979. 

J.  M.  Hovem,  "Finite  Amplitude  Effects  in  Marine  Sediments,"  Applied  Re¬ 
search  Laboratories  Technical  Paper  79-41  (ARL-TP-79-41 ) ,  Applied  Research 
Laboratories,  The  University  of  Texas  at  Austin,  June  1979. 

J.  M.  Hovem,  "Finite  Amplitude  Effects  in  Marine  Sediments,"  Applied  Re¬ 
search  Laboratories  Technical  Oral  Presentation  No.  79-7  (ARL-T0-79-7) , 

Applied  Research  Laboratories,  The  University  of  Texas  at  Austin,  September 
1979. 


56 


J.  M.  Hovem,  "Viscous  Attenuation  of  Sound  in  Suspensions  and  High  Porosity 
Marine  Sediments,"  Applied  Research  Laboratories  Technical  Paper  No.  79- 
62  (ARL-TP-79-62),  Applied  Research  Laboratories,  The  University  of  Texas 
at  Austin,  September  1979. 

J.  M.  Hovem,  "Finite  Amplitude  Effects  in  Marine  Sediments,"  presented 
at  the  Conference  on  Underwater  Applications  of  Nonlinear  Acoustics,  Uni¬ 
versity  of  Bath,  Bath,  England,  September  1979. 

J.  M.  Hovem,  "The  Nonlinearity  Parameter  of  Saturated  Marine  Sediments," 

J.  Acoust.  Soc.  Am.  66,  1463-1467  (1979). 

J.  M.  Hovem,  "Viscous  Attenuation  of  Sound  in  Suspensions  and  High  Porosity 
Marine  Sediments,"  J.  Acoust.  Soc.  Am.  62,  1559-1563  (1980). 

J.  M.  Hovem  and  G.  0.  Ingram,  "Viscous  Attenuation  of  Sound  in  Saturated 
Sand,"  Applied  Research  Laboratories  Technical  Paper  No.  79-35  (ARL-TP- 
79-35),  Applied  Research  Laboratories,  The  University  of  Texas  at  Austin, 
April  1979. 

J.  M.  Hovem  and  G.  D.  Ingram,  "Viscous  Attenuation  of  Sound  in  Saturated 
Sand,"  J.  Acoust.  Soc.  Am.  66^  1807-1812  (1979). 

A.  C.  Kibblewhite,  "Attenuation  of  Underwater  Sound  of  Low  Frequencies," 
Applied  Research  Laboratories  Technical  Report  No.  76-1  (ARL-TR-76-1), 
Applied  Research  Laboratories,  The  University  of  Texas  at  Austin,  December 
1975. 

D.  J.  Shirley,  "Final  Report  under  Contract  N00014-70-A-0166,  Task  0005," 
Applied  Research  Laboratories  Technical  Report  No.  72-6  (ARL-TR-72-6) , 
Applied  Research  Laboratories,  The  University  of  Texas  at  Austin,  January 
1972. 

D.  J.  Shirley,  "Interim  Technical  Description  of  the  ARL  Compressional 
Wave  In  Situ  Core  Profilometer,"  Applied  Research  Laboratories  Technical 
Memorandum  No.  74-9  (ARL-TM-74-9) ,  Applied  Research  Laboratories,  The  Uni¬ 
versity  of  Texas  at  Austin,  March  1974. 

D.  J.  Shirley,  "Calibration  Manual  for  ARL  Profilometer,"  Informal  Memoran¬ 
dum,  July  1974. 

D.  J.  Shirley,  "Fine  Structure  of  the  Sound  Speed  Profile  in  Ocean  Bottom 
Sediments  from  In  Situ  Measurements,"  presented  at  the  89th  Meeting  of 
the  Acoustical  Society  of  America,  Austin,  Texas,  8-11  April  1975. 

D.  J.  Shirley,  "Transducer  for  Generation  and  Detection  of  Shear  Waves," 

ARL  Invention  Disclosure,  September  1975. 

D.  J.  Shirley,  "Determination  of  the  Acoustic  Properties  of  Deep  Ocean 
Sediments  from  In  Situ  Profiles,"  presented  at  the  92nd  Meeting  of  the 
Acoustical  Society  of  America,  San  Diego,  California,  16-19  November  1976. 


D.  J.  Shirley,  "Method  for  Measuring  In  Situ  Acoustic  Impedance  of  Marine 
Sediments,"  Applied  Research  Laboratories  Technical  Paper  No.  77-19  (ARL- 
TP-77-19),  Applied  Research  Laboratories,  The  University  of  Texas  at  Austin, 
May  1977. 

0.  J.  Shirley,  "Acoustic  Impedance  Measuring  Device  for  Marine  Sediments," 
presented  at  the  93rd  Meeting  of  the  Acoustical  Society  of  America,  Univer¬ 
sity  Park,  Pennsylvania,  6-10  June  1977. 

D.  J.  Shirley,  "Laboratory  and  In  Situ  Sediment  Acoustics,"  Applied  Research 
Laboratories  Technical  Report  No.  77-46  (ARL-TR-77-46) ,  Applied  Research 
Laboratories,  The  University  of  Texas  at  Austin,  August  1977. 

D.  J.  Shirley,  "An  Improved  Shear  Wave  Transducer,"  Applied  Research  Labora¬ 
tories  Technical  Paper  No.  77-39  (ARL-TP-77-39) ,  Applied  Research  Labora¬ 
tories,  The  University  of  Texas  at  Austin,  December  1977. 

D.  J.  Shirley,  "Method  for  Measuring  In  Situ  Acoustic  Impedance  of  Marine 
Sediments,"  J.  Acoust.  Soc.  Am.  62_,  1028-1032  (1977). 

D.  J.  Shirley,  "An  Improved  Shear  Wave  Transducer,"  J.  Acoust.  Soc.  Am. 

63,  1643-1645  (1978). 

D.  J.  Shirley,  "The  ARL:UT  Subseafloor  Environmental  Simulator,"  Applied 
Research  Laboratories  Technical  Report  No.  79-53  (ARL-TR-79-53) ,  Applied 
Research  Laboratories,  The  University  of  Texas  at  Austin,  November  1979. 

D.  J.  Shirley,  "Acoustical  Properties  of  Sediments,"  Applied  Research  Lab¬ 
oratories  Technical  Report  No.  81-20  (ARL-TR-81-20) ,  Applied  Research  Lab¬ 
oratories,  The  University  of  Texas  at  Austin,  May  1981. 

D.  J.  Shirley  and  A.  L.  Anderson,  "Compressional  Wave  Profilometer  for 
Deep  Water  Measurements,"  Applied  Research  Laboratories  Technical  Report 
No.  74-51  (ARL-TR-74-51 ) ,  Applied  Research  Laboratories,  The  University 
of  Texas  at  Austin,  December  1974. 

D.  J.  Shirley,  and  A.  L.  Anderson,  "Studies  of  Sediment  Shear  Waves,  Acous¬ 
tical  Impedance,  and  Engineering  Properties,"  Applied  Research  Laboratories 
Technical  Report  No.  75-23  (ARL-TR-75-23),  Applied  Research  Laboratories, 

The  University  of  Texas  at  Austin,  May  1975. 

D.  J.  Shirley  and  A.  L.  Anderson,  "Acoustic  and  Engineering  Properties 
of  Sediments,"  Applied  Research  Laboratories  Technical  Report  No.  75-58 
(ARL-TR-75-58) ,  Applied  Research  Laboratories,  The  University  of  Texas 
at  Austin,  October  1975. 

D.  J.  Shirley  and  A.  L.  Anderson,  "In  Situ  Measurement  of  Marine  Sediment 
Acoustical  Properties  During  Coring  in  Deep  Water,"  Applied  Research  Labora¬ 
tories  Technical  Paper  No.  75-3  (ARL-TP-75-3) ,  Applied  Research  Laboratories, 
The  University  of  Texas  at  Austin,  October  1975. 


58 


D.  J.  Shirley  and  A.  L.  Anderson,  "In  Situ  Measurement  of  Marine  Sediment 
Acoustical  Properties  During  Coring  in  Deep  Water,"  IEEE  Transactions  on 
Geoscience  Electronics,  Vol .  GE-13,  No.  4,  pp.  163-169,  October  1975. 

D.  J.  Shirley  and  A.  L.  Anderson,  "Experimental  Investigation  of  Shear 
Waves  in  Laboratory  Sediments,"  presented  at  the  90th  Meeting  of  the  Acous¬ 
tical  Society  of  America,  San  Francisco,  California,  3-7  November  1975. 

D.  J.  Shirley  and  A.  L.  Anderson,  "Shear  Waves  in  Unconsolidated  Sediments," 
presented  at  the  92nd  Meeting  of  the  Acoustical  Society  of  America,  San 
Diego,  California,  16-19  November  1976. 

D.  J.  Shirley,  A.  L.  Anderson,  and  L.  D.  Hampton,  "In  Situ  Measurement 
of  Sediment  Sound  Speed  During  Coring,"  Applied  Research  Laboratories  Tech¬ 
nical  Report  No.  73-1  (ARL-TR-73-1 ) ,  Applied  Research  Laboratories,  The 
University  of  Texas  at  Austin,  March  1973. 

D.  J.  Shirley,  A.  L.  Anderson,  and  L.  D.  Hampton,  "Measurement  of  In  Situ 
Speed  During  Sediment  Coring,"  OCEAN  *73,  Record  of  the  International  Con¬ 
ference  on  Engineering  in  the  Ocean  Environment,  Seattle,  Washington,  25- 
28  September  1973. 

D.  J.  Shirley,  A.  Bedford,  and  S.  K.  Mitchell,  "Wave  Propagation  in  a  Satu¬ 
rated  Model  Sediment  with  Varied  Liquid  Properties,"  Submitted  to  J.  Acoust. 
Soc.  Am.  (1982). 

D.  J.  Shirley  and  D.  W.  Bell,  "Acoustics  of  In  Situ  and  Laboratory  Sedi¬ 
ments,"  Applied  Research  Laboratories  Technical  Report  No.  78-36  (ARL-TR- 
78-36),  Applied  Research  Laboratories,  The  University  of  Texas  at  Austin, 
August  1978. 

D.  J.  Shirley  and  D.  W.  Bell,  "Temperature  Variation  of  the  Acoustic  Proper¬ 
ties  of  Laboratory  Sediments,"  presented  at  the  98th  Meeting  of  the  Acous¬ 
tical  Society  of  America,  Salt  Lake  City,  Utah,  November  1979. 

D.  J.  Shirley,  D.  W.  Bell  and  J.  M.  Hovem,  "Laboratory  and  Field  Studies 
of  Sediment  Acoustics,"  Applied  Research  Laboratories  Technical  Report 
No.  79-26  (ARL-TR-79-26) ,  Applied  Research  Laboratories,  The  University 
of  Texas  at  Austin,  June  1979. 

D.  J.  Shirley  and  L.  D.  Hampton,  "Acoustic  Velocity  Profilometer  for  Sedi¬ 
ment  Cores,"  OCEAN  "72,  Record  of  the  International  Conference  on  Engineer¬ 
ing  in  the  Ocean  Environment,  Newport,  Rhode  Island,  13-15  September  1972. 

D.  J.  Shirley  and  L.  D.  Hampton,  "Determination  of  Sound  Speed  in  Cored 
Sediments,"  Applied  Research  Laboratories  Technical  Report  No.  72-44  (ARL- 
TR-72-44),  Applied  Research  Laboratories,  The  University  of  Texas  at  Austin, 
December  1972. 

D.  J.  Shirley  and  L.  D.  Hampton,  "Acoustic  Velocimeter  for  Ocean  Bottoms 
Coring  Apparatus,"  ARL  Invention  Disclosure,  November  1973. 


59 


D.  J.  Shirley  and  L.  0.  Hampton,  "Shear  Wave  Measurements  in  Laboratory 
Sediments,"  Applied  Research  Laboratories  Technical  Paper  No.  77-31  ( ARL- 
TP-77-31),  Applied  Research  Laboratories,  The  University  of  Texas  at  Austin, 
August  1977. 

D.  J.  Shirley  and  L.  D.  Hampton,  "Shear  Wave  Measurements  in  Laboratory 
Sediments,"  J.  Acoust.  Soc.  Am.  6^.  607-613  (1978). 

D.  J.  Shirley,  J.  M.  Hovem,  G.  D.  Ingram,  and  D.  W.  Bell,  "Sediment  Acous¬ 
tics,"  Applied  Research  Laboratories  Technical  Report  No.  80-17  (ARL-TR- 
80-17),  Applied  Research  Laboratories,  The  University  of  Texas  at  Austin, 
April  1980. 

K.  H.  Stokoe  (UT  Austin),  D.  G.  Anderson  (Fugro,  Inc.),  E.  J.  Arnold  (UT 
Austin),  R.  J.  Hoar  (UT  Austin),  and  0.  J.  Shirley,  "Development  of  a  Bottom 
Hole  Device  for  Offshore  Shear  Wave  Velocity  Measurement,"  Proceedings 
of  the  10th  Annual  Offshore  Technology  Conference,  Houston,  Texas,  8-11 
May  1978. 


B.  E.  Tucholke  (Lamont-Doherty  Goelogical  Observatory,  Columbia  University) 
and  D.  J.  Shirley,  "Comparison  of  Laboratory  and  In  Situ  Congressional  - 
Wave  Velocity  Measurements  on  Sediment  Cores  from  the  Western  North  Atlantic," 
Journal  of  Geophysical  Research,  84,  687-695  (1979). 


60 


2  September  1982 


DISTRIBUTION  LIST  FOR 
ARL-TR-82-47 

UNDER  CONTRACT  N00014-76-C-0117 
UNCLASSIFIED 


Copy  No. 


1 

2 

3 

4 

5 

6 


7 

8 
9 

10 


11 

12 

13 


14 

15 


Commanding  Officer 

Naval  Ocean  Research  and  Development  Activity 
NSTL  Station,  MS  39529 
Attn:  Code  110 

W.  Kuperman  (Code  320) 

Code  340 

R.  Gardner  (Code  520) 

E.  Chaika  (Code  530) 

H.  Eppert  (Code  300) 

Commanding  Officer 

Office  of  Naval  Research 

Arlington,  VA  22217 

Attn:  CAPT  A.  Gilmore  (Code  200) 

J.  McKisic  (Code  325A) 

R.  Obrochta  (Code  464) 

CDR  M.  McDonald 

Commander 

Naval  Electronic  Systems  Command 
Department  of  the  Navy 
Washington,  DC  20360 
Attn:  J.  Reeves  (Code  PME  124-30) 

CDR  C.  Spikes  (Code  PME  124-60) 

LCDR  S.  Hollis  (Code  6121) 

Commander 

Naval  Sea  Systems  Command 
Department  of  the  Navy 
Washington,  DC  20362 
Attn:  R.  Farwell  (Code  63R) 

Chief  of  Naval  Material 
Office  of  Naval  Technology 
Department  of  the  Navy 
Washington,  DC  20360 


61 


Distribution  List  for  ARL-TR-82-47  under  Contract  N00014-76-C-0117 
(Cont'd) 

Copy  No. 

Commanding  Officer 

Naval  Oceanographic  Office 

NSTL  Station,  Bay  St.  Louis,  MS  39522 

16  Attn:  W.  Jobst 

17  M.  G.  Lewis 

Commander 

Naval  Ocean  Systems  Center 
Department  of  the  Navy 
San  Diego,  CA  92132 
Attn:  Library 

E.  L.  Hamilton 
H.  P.  Bucker 

Director 

Naval  Research  Laboratory 
Department  of  the  Navy 
Washington,  DC  20375 
Attn :  B .  Adams 

R.  Moseley 
Code  2627 

Chief  of  Naval  Operations 
Department  of  the  Navy 
Washington,  DC  20350 

27  Attn:  CAPT  E.  Young  (OP-9521) 

Commander 

Naval  Air  Development  Center 
Department  of  the  Navy 
Warminster,  PA  18974 

28  Attn:  C.  L.  Bartberger 

Commander 

New  London  Laboratory 
Naval  Underwater  Systems  Center 
Department  of  the  Navy 
New  London,  CT  06320 

29  Attn:  B.  Cole  (Code  24) 

Commander 

Naval  Coastal  Systems  Center 
Panama  City,  FL  32401 

30  Attn:  E.  G.  McLeroy,  Jr. 

31  B.  Tolbert 


62 


Distribution  List  for  ARL-TR-82-47  under  Contract  N00014-76-C-0117 
(Cont'd) 

Copy  No. 

Commanding  Officer 

Naval  Facilities  Engineering  Command 
Department  of  the  Navy 
200  Stovall  Street 
Alexandria,  VA  22332 

32  Attn:  R.  A.  Peloquin,  Code  0320 
Commanding  Officer 

U.S.  Naval  Civil  Engineering  Laboratory 
Port  Hueneme,  CA  93043 

33  Attn:  N.  Albertson 

Superintendent 

Naval  Postgraduate  School 

Monterey,  CA  93940 

34  Attn:  H.  Medwin,  Physics  Dept. 

35  Library 

Woods  Hole  Oceanographic  Institution 
Woods  Hole,  MA  02543 
Zo  ..  Attn:  C.  Hollister 

37  B.  Tucholke 

38  E.  Hayes 

Hawaii  Institute  of  Geophysics 
The  University  of  Hawaii 
2525  Correa  Road 
Honolulu,  HI  96822 

39  Attn:  G.  Sutton 

40  M.  Manghnani 

The  Scripps  Institution  of  Oceanography 
The  University  of  California/San  Diego 
San  Diego,  CA  92151 

41  Attn:  P.  Lonsdale 

42  R.  Tyce 

43  F.  Fisher 

Department  of  Geological  Oceanography 
Texas  A&M  University 
College  Station,  TX  77840 

44  Attn:  W.  R.  Bryant 

Geophysics  Laboratory 
Marine  Science  Institute 
The  University  of  Texas  at  Austin 
700  The  Strand 
Galveston,  TX  77550 
Attn:  E.  W.  Behrens 


45 


Distribution  List  for  ARL-TR-82-47  under  Contract  N00014-76-C-0117 
(Cont'd) 

Copy  No. 

The  Catholic  University  of  America 
6220  Michigan  Avenue,  NE 
Washington,  DC  20017 

46  Attn:  H.  M.  Uberall 

Lamont-Doherty  Geological  Observatory 
Palisades,  NY  10964 

47  Attn:  G.  Bryan 

48  W.  J.  Ludwig 

49  R.  D.  Stoll 

Department  of  Civil  and  Ocean  Engineering 
The  University  of  Rhode  Island 
Kingston,  RI  02881 

50  Attn:  A.  J.  Silva 

University  College  of  North  Wales 
Marine  Science  Laboratories 
Menai  Bridge 
Anglesey,  NORTH  WALES 

51  Attn:  D.  Taylor  Smith 

52  P.  Schultheiss 

Department  of  Physics 
The  University  of  Auckland 
Auckland,  NEW  ZEALAND 

53  Attn:  A.  Kibblewhite 

Department  of  Civil  Engineering 
The  University  of  Texas  at  Austin 
Austin,  TX  78712 

54  Attn:  K.  Stokoe 

Southwest  Research  Institute 
P.  0.  Drawer  28510 
San  Antonio,  TX  78284 

55  Attn:  T.  Owen 

56  D.  J.  Shirley 

Western  Electric  Company 
Department  7426 
P.  0.  Box  20046 
Greensboro,  NC  27420 
Attn:  T.  F.  Clark 


57 


Distribution 

(Cont'd) 

Copy  No. 


58 

59 

60 


61 

62  -  73 


74 

75 

76 


77 


List  for  ARL-TR-82-47  under  Contract  N00014-76 


Electronics  Research  Laboratory 

The  University  of  Trondheim 

The  Norwegian  Institute  of  Technology 

O.  S.  Bragstad  Plass  6 
N-7034  Trondheim  -  NTH,  NORWAY 
Attn:  J.  Hovem 

Applied  Physics  Laboratory 
The  Johns  Hopkins  University 
Johns  Hopkins  Road 
Laurel,  MD  20812 
Attn:  J.  Lombardo 

Applied  Research  Laboratory 
The  Pennsylvania  State  University 

P.  0.  Box  30 

State  College,  PA  16801 

Department  of  Geology 
The  University  of  Texas  at  Austin 
Austin,  TX  78712 
Attn:  M.  Backus 

Commanding  Officer  and  Director 
Defense  Technical  Information  Center 
Cameron  Station,  Building  5 
5010  Duke  Street 
Alexandria,  VA  22314 

Director 

SACLANT  ASW  Research  Centre 
La  Spezia,  ITALY 
Attn:  R.  Goodman 

Defence  Research  Establishment  Atlantic 
9  Grove  Street 
Dartmouth,  N.S.,  CANADA 
Attn:  Library 

Office  of  Naval  Research 
Branch  Office  Chicago 
Department  of  the  Navy 
536  S.  Clark,  Room  286 
Chicago,  IL  60605 

Office  of  Naval  Research 
Resident  Representative 
Room  582,  Federal  Building 
Austin,  TX  78701 


-C-0117 


65 


Distribution  List  for  ARL-TR-82-47  under  Contract  N00014-76-C-0117 
(Cont'd) 

Copy  No. 

78  Environmental  Sciences  Division,  ARL : UT 

79  David  T.  Blackstock,  ARL:UT 

80  Harlan  G.  Frey,  ARL:UT 

81  Stephen  K.  Mitchell,  ARLtUT 

82  Thomas  G.  Muir,  ARL:UT 

83  Library,  ARL:UT 

84  -  90  Reserve,  ARL:UT 


